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We theoretically model the electronic dynamics of a coupled quantum dot pair in a static elec¬ 
tric field. We then investigate the possibility of polarization-entangled photon emission from the 
radiative cascade of the molecular biexciton state. Through numerical simulations, we analyze the 
dependence of entanglement fidelity on temperature and electric field, as well as tunnel coupling. 
We establish a regime of direct-indirect exciton detunings for which coupled quantum dots are su¬ 
perior to single dots for entangled photon generation, yielding near-unit fidelity over a larger range 
of exchange splittings. 


I. INTRODUCTION 

Quantum mechanics allows for the unique possibility of 
instantaneous non-local correlation, or entanglement, be¬ 
tween particles. Several technologies have emerged which 
use quantum entanglement as a fundamental resource, in¬ 
cluding quantum computation 1 ^ and cryptography,^^ 
as well as entangled-photon microscopy. 11 Quantum in¬ 
formation technologies can be implemented by using the 
orthogonal polarization states of a photon to encode bi¬ 
nary information in a “dual-rail” qubit scheme,^ which 
allows for simple manipulation and distribution of pho¬ 
tonic qubits using linear optical elements,^ 3 or using spin 
states of bound charges P^l However, a source of entan¬ 
glement between photons or distant spins is needed to 
implement two-qubit gates or two-photon quantum cryp¬ 
tography protocols. 

Entangled photon pairs are routinely produced by non¬ 
linear optical processes such as spontaneous parametric 
downconversion (SPDC), in which one photon interacts 
with a bulk crystal to produce two lower-energy photons 
with opposite polarizations. However, due to the prob¬ 
abilistic nature of such an interaction, SPDC-based en¬ 
tangled photon sources can always produce zero, one, or 
multiple photon pairs from each incident pulse, which 
can introduce errors in entanglement-based quantum 
information! 1 ^ Semiconductor quantum dots (QDs) are 
promising candidates for on-denrand sources of enta ngled 
photon pairs, as well as spin-photon entanglement 
Two-photon emission can be achieved by optically excit¬ 
ing the neutral biexciton state, consisting of two bound 
excitons. Radiative decay from this state can proceed 
along one of two pathways corresponding to the two op¬ 
tically active single-exciton spin states. Coherence be¬ 
tween these two oppositely-polarized (but otherwise in¬ 
distinguishable) emission processes results in a pair of 
polarization-entangled photons. In single QDs, however, 
the degeneracy of the intermediate exciton spin states is 
often removed by the anisotropic electron-hole exchange 
splitting, which is due to asymmetry in the self-assembled 
QD system PSHU] jf this splitting is larger than the opti¬ 
cal transition linewidths, the emitted photons are dis¬ 
tinguishable by their energy and exhibit correlated—but 
not entangled -polarizations? 21 Several techniques have 


been implemented to avoid this problem, including post¬ 
selection via spectral filtering 22 ! and minimization of the 
exchange splitting by postgr owt h thermal annealingpU or 
application of external fieldsP® 5 Alternative QD growth 
methods have also been demonstrated to exhibit a sys¬ 
tematically smaller exchange splitting, including inter¬ 
face fluctuation QDsP and QDs grown on the (111) sur¬ 
face of GaAs! 2 '-^ Aside from neutral biexciton decay, 
negatively-charged states in QDs can also be used to cre¬ 
ate a chain of entangled photons in a one-dimensional 
cluster stated This procedure can be generalized to two 
dimensions using two vertically-stacked coupled quantum 
dots (CQDs)P^Here we discuss a scheme in which entan¬ 
gled photon pairs are obtained straightforwardly through 
a cascaded biexciton-exciton decay. 

Our scheme for on-denrand entangled photon pairs uses 
CQDs grown in an electric field-effect structure, such as 
a Schottky diode!^2ll The CQD structure permits two 
types of neutral excitons: direct excitons, with the elec¬ 
tron and hole in the same dot, and indirect excitons, 
with the electron and hole separated in different dots. 
Because of their spatial charge separation, the energies 
of indirect states can easily be tuned by applying an 
electric field. 35 3 ' In addition, spatial charge separation 
greatly reduces both the short-range isotropic and long- 
range anis otropi c exchange splitting of indirect exciton 
spin statesP 3 38 -! This property can potentially be used 
for generation of entangled photon pairs by preparing 
the CQD in an indirect biexciton spin singlet state, with 
one direct (single-dot) exciton and one indirect exciton. 
In an effort to determine the utility of this proposal, we 
theoretically calculate the fidelity of entangled photon 
generation in a CQD system. 

II. THEORETICAL MODEL 

We consider an asymmetric CQD system of two verti¬ 
cally stacked, tunnel-coupled, self-assembled InAs-GaAs 
QDs (referred to here as the bottom and top dots) with 
respective heights hg and h t and center-to-center sepa¬ 
ration d , grown in a diode structure to allow application 
of an electric field F along the growth direction. Depend¬ 
ing on the dot asymmetry and applied field direction, ei¬ 
ther electron or hole levels can be tuned into resonance, 
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FIG. 1: (a) Energy level diagram of the molecular biexciton radiative cascade, depicting polarizations of all allowed recom¬ 
bination pathways, as well as spin-flip and phonon-assisted hole tunneling channels. Thick vertical arrows indicate direct 
exciton recombination, while thin vertical arrows indicate a slower interdot recombination. Excitation and detection pathways 
are indicated by green shading, (b) Simulated electric field-dispersed photoluminescence spectrum, mapping out each of the 
relevant transition energies in the vicinity of the exciton and biexciton hole tunneling resonances. Single-exciton hole tunneling 
resonance occurs at A = 0 meV, while biexciton hole tunneling resonances occur at A = —1.9 and 27.2 meV. 


causing electrons or holes to tunnel coherently.®! For the 
remainder of this work and without loss of generality we 
will assume a CQD diode structure is chosen that pro¬ 
motes hole tunneling, with any excited electrons confined 
to the bottom dot. 

We restrict the CQD Hilbert space to that of the crys¬ 
tal ground state | g), direct exciton |X), indirect exci¬ 
ton \iX), and molecular biexciton \iXX) (consisting of 
one direct and one indirect exciton). Each of the single- 
exciton charge states can exist in one of four spin states: 
two optically active bright states | H/V) = llff ± til), 
and two optically inactive dark states \H d /Vd) = 
|tfl ± ID) (single (double) arrows denote spin-1/2 elec¬ 
trons (spin-3/2 heavy holes)). The molecular biexci¬ 
ton can exist in a singlet | S) = |f| — It) Ill'll — Jill) or 
one of three triplet spin states |T + ) = |f| — It) lltH), 

\T 0 ) = |tl - It) \U + U), or |T_) = |tl - It) |DD)- 

While a CQD exhibits many optical transitions, we 
focus on the molecular biexciton spin singlet state 
| iXX,S), which can produce correlated photon pairs by 
sequential recombination to the ground state due to the 
optical selection rules. As depicted in Fig. [lja), this 
state is optically coupled to the ground state through the 
bright spin states of the direct and indirect exciton, form¬ 
ing two pairs of oppositely-polarized decay pathways. 

Resonant excitation can selectively prepare the system 
in the molecular biexciton state. In contrast to biex- 
citons in single QDs, the molecular biexciton can be 
resonantly excited by stepwise two-laser driving along 
one of two pathways, via the direct or indirect exciton 
state. Since the indirect exciton pathway features a much 
smaller exchange splitting, the preferred two-laser exci¬ 
tation scheme is to drive each transition of the direct ex¬ 


citon pathway and detect the emission along the indirect 
exciton pathway, as depicted in Fig. [lja). The relevant 
transitions are identified in the simulated electric field- 
dispersed photoluminescence spectrum of Fig. Bb). The 
singlet spin state can be spectrally selected near a biex¬ 
citon tunneling resonance, with triplet states excluded 
from the spin-conserving tunnel coupling, or by utiliz¬ 
ing two-photon absorption, due to the selection rules of 
simultaneous transitions.®! 

We consider the case of pulsed excitation, which al¬ 
lows for on-demand preparation of the molecular biex¬ 
citon state and subsequent emission of a single entan¬ 
gled photon pair. For two-laser driving, one must fur¬ 
ther specify the relative timing of the two pulses: they 
may be applied sequentially, or concurrently. Bensky et 
alpU examined each of these pulsed excitation schemes in 
the case of a single QD and found that both can gener¬ 
ate the biexciton state with more than 99% probability, 
though the concurrent pulse scheme may be less reliable 
due to interference between the pulses. Pulsed excitation 
of the molecular biexciton has not been examined exper¬ 
imentally, though continuous-wave excitation has been 
demonstrated.®! 

Photon emission events can be described within the 
CQD and photon subspaces using transition operators 
< Ti h/v,d/i ( &2H/v,d/i ) and photon creation operators 
a ' H / v > which act on the system to produce the first 
(second) ff/H-polarized photon by direct/indirect re¬ 
combination of the biexciton (exciton) state (see Ta¬ 
ble I for details). Assuming that the system has been 
prepared in the molecular biexciton spin singlet state 
|\P 0 ) = \iXX, S) and that recombination occurs via the 
indirect exciton pathway without nonradiative transi- 
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tions, the first recombination produces the state 


’hi) = + flydi v,d) I’Fo) 

= ^=(\H) 1 \iX,H) + \V) 1 \iX,V)). 


(1) 


During the time r before the second photon is emit¬ 
ted, the system undergoes coherent evolution, acquiring 
a phase difference due to the exchange interaction: 


|H/ 1 (r))=e- i ^/ fi |^ 1 (0)) 

= \iX, H) + e iSlT / h |F) 1 1 iX , V)), 

( 2 ) 

to within a global phase factor, where Si is the exchange 
splitting between indirect exciton bright spin states. Fi¬ 
nally, the second recombination results in 


1^2) = + a\f(t2v,i) I’M 7 ")) 

= ^=(\H) 1 \H) 2 + e iS ^ h \V) 1 \V) 2 )\g) 


( 3 ) 


= |tf 


2 p, 


1 9 ) 


with the two photons in a pure polarization-entangled 
state | v I , 2p). However, since recombination can occur via 
multiple pathways, the result is a mixed state, described 
by the two-photon density matrix p> 2 P- 

In the absence of anisotropic electron-hole exchange 
splitting, each pair of oppositely-polarized decay path¬ 
ways are energetically indistinguishable, resulting in the 
maximally-entangled Bell state |’F + ) = (| H) 1 \H) 2 + 
|V)i \V) ^)/-\/ 2 for the polarizations of the emitted pho¬ 
ton pair. 42 Under generic conditions, we quantify the 
quality of entanglement by the fidelity of the two-photon 
polarization state to the target Bell state. This fidelity 
is determined by the projection 


F + = (T+l /5 2 p |’F + ) 

= 2 ^P\ hh )( hh \ + P\VV)(VV\) + R- e (P\HH)(VV\)- 

(4) 

A fidelity higher than 0.5 indicates polarization entan¬ 
glement, while a lower fidelity is acheivable by classical 
correlation alone. Uncorrelated photons, such as back¬ 
ground light not o rigina ting from the CQD, give a base¬ 
line fidelity of 0.25 .®® 

The CQD system interacts with the finite-temperature 
phonon bath of the surrounding crystal lattice, as well as 
the vacuum electromagnetic field via spontaneous emis¬ 
sion. We assume Markovian behavior, such that the envi¬ 
ronment always remains at equilibrium (emitted particles 
cannot be reabsorbed). This assumption requires that 
any system-environment correlations decay on a much 
faster time scale than the system dynamics, which could 
introduce some inaccuracies at the lowest temperatures 


TABLE I: List of incoherent processes included in the master 
equation, with corresponding transition rates and operators. 
Exciton spin states are indexed by m = H, V, H d , V d . 


Process 

Rate 

Operator 

Recombination 

Tx 

&m/v,D = \iX,H/V) {iXX,S\ 



*2H/V,D = 1 9) ( iX,H/V\ 


r ix 

&1H/VJ = \X,H/V) (iXX,S\ 



*2h/v,i = Iff) (iX,H/V | 

Phonon-assisted tunneling 

r abs 

| A', m) (iX, m\ 


r 

em 

| iX, m) {X, m\ 

Phonon-assisted spin-flip 

1 Sf 

| X,H d /V d ) {X,H/V |, H.c. 


piX 

1 sf 

| iX,H d /V d ) {iX,H/V |, H.c. 

Pure dephasing 

7 p 

E m |A,m) (X, m\ 


Y P X 

E m | iX,m) (iX,m\ 


simulated (T = 1 K) due to long-lived reservoir correla¬ 
tions. Using the Rotating Wave Approximation to aver¬ 
age over quickly-oscillating terms, we can then describe 
the dissipative dynamics of this open quantum system 
using the Lindblad master equation^ for the reduced 
density matrix p of the system: 


= ±[H,p] + fajal - \a\aj \pa\a^ 

i x 7 

( 5 ) 

We assume that an initial state of p(0) = 

| iXX,S) (iXX,S\ is prepared by pulsed excitation, and 
do not explicitly include interaction with the optical 
pulse. By projecting onto the basis states, Eq. © can 
be written as (i\p\j) = pij = Yjki M ijMPkh where the 
elements of the time-dependence tensor M 13k i contain 
the various coherent phase evolution, transition and de¬ 
cay rates. Vectorizing the elements of p results in the 
matrix differential equation p = Mp, leading to the so¬ 
lution p(t) = e Mt p( 0). In terms of the superoperator 
J 2 fp = d t p defined by Eq. ([ 5 ]), this solution can also be 
written as p(t) = e^ T p( 0). Details of the numerical sim¬ 
ulation are given in the appendix, including calculation 
of single-particle wavefunctions, Hamiltonian matrix ele¬ 
ments, and transition rates. 

For a given decay pathway, with the first (second) pho¬ 
ton generated by a((3) = D/I recombination after a time 
delay r, the density matrix element (ij \ fop (t) | kl) in 
the two-photon linear polarization basis is related to the 
CQD dynamics by the two-time correlation functiorP^ 47 ^ 


ffijh, a p( T ) = ( a ii,a(°) a 2 j,p( T ) & 2 l,p(r)a lkta ( 0 )) 
= Tr e^ T {ai kta p(0)a\ i a ) a\ j p 


- /At 


( 6 ) 


where the last line follows from the quantum regression 
theorem.® These polarization correlations can be mea¬ 
sured experimentally using time-correlated single-photon 
counting in a Hanbury Brown-Twiss setup,® where 
each detection pathway contains polarization optics 
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and a monochromator to spectrally select the relevant 
transitions and perform polarization cross-correlation 
measurements.^ 2 42 Since recombination can occur via 
several pathways with a random time delay, the elements 
of p 2 p are calculated by using a time-averaged statistical 
mixture 

p\ij ){ ki\=A J2 [ dTP ^P2p(T)g^l lafi (T), (7) 

a,/3=D,I J ° 


where Td is the detection time window, A is a normaliza¬ 
tion constant set to enforce the condition Tr(p 2 p) = 1, 
P\D(i) = r_Y(i;v)/(rx + IYy) is the probability of direct 
(indirect) recombination of the first exciton, P 2 D(i){t ) = 
^x(iX)nx{iX)( T ) is the probability per unit time of the 
second direct (indirect) recombination, and rixux)( T ) is 
the population of bright states |X(?X),±1) at time r. 
Transition rates for direct and indirect exciton recombi¬ 
nation are denoted by Ty and Tix, respectively. 

The critical dynamics determining photon entangle¬ 
ment occurs during the time r between photon emission 
events, when the dynamics is limited to the single-exciton 
subspace. Truncating to the bright spin states, the sim¬ 
plified exciton Hamiltonian can be written 


H = 



1 iX,V) 

IX, H) 

\X,V) 

0 

th 

0 

-Si /2 

0 

th 

0 

A + S D / 2 

0 

th 

0 

A - S d /2 


( 8 ) 

where A = ed(F — F 0 ) is the exciton detuning due to the 
applied electric field, Fo is the field value at the |X) — \iX) 
anticrossing, and th is the resonant hole tunnel coupling 
including Coulomb correction. The magnitude of the ex¬ 
change splitting S D /i is a critical factor for entangled 
photon emission, as a splitting larger than the radiative 
linewidth of the corresponding exciton is expected to ren¬ 
der the two decay pathways distinguishable and prevent 
entanglement. The bright-state exchange splitting is ap¬ 
proximated using the dipole interaction term of the long- 
range multipole expansion®! 


S D(i) « — 


Prd 6 




|r — r 

l,h* 


'13 


( 9 ) 


X V’f} ( r ) % l ) B(T)i r )V , B(T)(C)’ 


where e is the average permittivity of InAs and GaAs, 
1 is the 3x3 identity matrix, n is a unit vector in the 
direction of (r — r'), and p ax = e {u% \ r is the inter¬ 
band transition dipole moment. Since Sp varies widely 
in self-assembled QDs in the range 0—100 peV depending 
on detailed growth conditions, we leave it as a variable 
parameter in the model and numerically evaluate Eq. 
to obtain Sj/Sd (see the appendix for details). 

Incoherent transitions, described by operators <7,; and 
occuring at rates Tj, include optical recombination, 
phonon-assisted hole tunneling, spin-flip, and pure de¬ 
phasing (see Table |I|. These transitions are incoherent in 


the sense that they occur at random times, because we 
trace over the photonic and phononic environment de¬ 
grees of freedom to obtain the reduced density matrix of 
the CQD system. We assume a direct exciton recombina¬ 
tion rate Ty = 1 ns -1 similar to experimentally observed 
rates in InAs-GaAs QDs. Indirect recombination rates 
Tix are found to be slower than direct recombination by 
a factor of ( Mbt/Mbb) 2 ~ 100 — 1000 due to a smaller 
electron-hole wavefunction overlap Mij = for 

the interdot state!® 

Holes undergo incoherent spin-conserving tunneling 
transitions via emission or absorption of a phonon, caus¬ 
ing transitions between direct and indirect exciton states 
at rates given by Fermi’s golden rule as 

2 7 r 

T abs = —J(A)n B (A,T) (10) 

2tt 

r em = -JT J(A)(n B {A, T) + 1), (11) 


where J{E) = | {X\ H e i~ ph \iX) \ 2 6(Eg’ - E) is 

the phonon spectral density of the interdot transition, 
Hb(E,T) = \/{e E P BT — 1) is the Bose distribution, giv¬ 
ing the population of phonon modes at energy E and tem¬ 
perature T, q denotes the phonon wave vector, Fl e i_ p h 
is the electron-phonon interaction Hamiltonian, and we 
use a linear dispersion relation E^ = Hcba |<?| for acous¬ 
tic phonons. We assume here that A = Ex — Eix > 0, 
so that tunneling to the direct exciton is partially sup¬ 
pressed, requiring absorption of a phonon. Since the de¬ 
tuning in the region of interest is much less than the 
LO phonon energy in GaAs, we consider only acoustic 
phonon-mediated transitions and neglect coupling to LO 
phonons. Including only interaction via the acoustic de¬ 
formation potential, which is usually found to dominate 
for electronic transitions in quantum dots, leads to the 
expression 51 


J{ A) = 


A a\ 


I67r 3 pc 2 


d 3 q 


LA 


^r^lBp^hT 


2 

S(E q ~ A), 


Q2) 

where ay is the valence-band deformation potential, p is 
the mass density of the crystal, c^a is the speed of the 
LA phonon mode in GaAs, and ipha is the single-particle 
wavefunction for a hole localized in dot a = B,T. We 
note that multiphonon and LO-phonon tunneling transi¬ 
tions are not included, which may become important at 
higher temperatures and exciton detunings than consid¬ 
ered here. 

Spin-flip transitions are also known to occur be¬ 
tween exciton spin states, via phonon and spin-orbit 
coupling!®!® j n particular, spin relaxation occurs in 
direct (indirect) excitons primarily via phonon-assisted 
transitions between bright and dark spin states at a 
rate Y^ tX> . Since our two-band effective mass model 
does not account for spin-orbit coupling, we include 
spin relaxation phenomenologically. We include only 
transitions between bright and dark spin states, since 
bright-bright and dark-dark transitions are at least an 
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FIG. 2: (a) Dependence of entanglement fidelity on direct 
exciton bright state splitting at T = 20 K, th = 0.2 peV, and 
A = 7 meV (green, top curve) and A = 3 meV (red, bottom 
curve), compared to the analogous SQD case, (b) Maximum 
fidelity and (c) width of F + (Sn) curve as a function of exciton 
detuning and temperature. 


order of magnitude slower due to their smaller tran¬ 
sition energies.^! We use the temperature-dependence 
of the spin-flip rate measured by Fras et alP^, scaled 
to fit the behavior measured by Hudson et alPS 
r ff(T) = (0.27 ns -1 ) + (0.29 /rs _ 1 K _ 2 )T 2 , correspond¬ 
ing to a thermally-activated two-phonon process. In the 
low-temperature regime where single-phonon transitions 
dominate and for small bright-dark splittings 
rf f (lX) is given by Fermi’s golden rule as in Eq’s (10) 
and {IlJ), with J( A° d {1) ) oc (Af d (/) ) 3 and n(Af d (/) ,T) « 
kjjT/A bd r> P-S Since the bright-dark splitting scales with 
electron-hole wavefunction overlap, we predict the rela¬ 
tionship = ( Mbt/Mbb ) 2 r between direct and in¬ 
direct spin-flip rates. 


Phonon-induced pure dephasing can have a substantial 
effect on interband coherences, most notably an increase 
of transition linewidth at elevated temperatures.^ Since 
phonons carry zero spin however, they couple equally 
to each electronic spin state. While, in principle, a 
slight spin-dependence of the phonon coupling could lead 
to a pure cross-dephasing 7 *J of the indirect exciton 
coherence,® we neglect this possibility for the remain¬ 
der of this work on the grounds that Moody et alP^ 
have measured the cross-dephasing rate to be negligible 
in the case of a single QD. As a result, the relative phase 
between exciton spin states—and the resulting photon 
entanglement—is largely unaffected by phonon-induced 
pure dephasing processes. 


III. ENTANGLEMENT 

The results of the numerical simulation are summa¬ 
rized in Fig. [2] The fidelity approximately exhibits a 
Lorentzian dependence on the direct exciton exchange 
splitting Sd, si milar to the behavior reported in the case 
of a single QD.EsHzHoHsl To characterize this behavior, 
we record the maximum fidelity F+ ax at Sd = 0 and the 
fidelity width AFg, defined as the largest value of \Sd\ 
with F + > .9. By plotting these fidelity characteristics 
as a function of the experimentally tunable exciton de¬ 
tuning A and temperature T, we obtain a map of possible 
behaviors observable in a single CQD pair. In general, 
we observe that the quality of entanglement depends pri¬ 
marily on the total dephasing rate 

7f = r iX + r abs + r iX (13) 

describing decay of the indirect exciton bright state co¬ 
herence (iX, +1| p | iX, —1), as well as the effective decay 
rate F e ff of the indirect exciton population (iX\ p\iX), 
determined by an exponential fit of the bright state den¬ 
sity matrix dynamics. For low temperatures (T < 40 K) 
and large detunings (A > 3 meV), entanglement fidelity 
can reach values above 0.99 for small exchange splittings, 
with a width A F~^ ss 1.5 — 2.7 peV. Note that the CQD 
system exhibits this high-fidelity behavior with Sd val¬ 
ues up to a factor of 5 greater than in the single QD 
case (A Fq S q D ss 0.4 — 0.6 peV). In this region, the 
decoherence is radiatively-limited (Y d x < 2 Fix)- allow¬ 
ing coherent evolution between photon emission events. 
With increasing lattice temperature, the overall fidelity 
decreases due to phonon-assisted tunneling and spin-flip. 
At detuning values below 3 meV, the large dephasing rate 
('J d x > 100 T ix ) prevents photon entanglement even in 
the absence of anisotropic exchange splitting, giving a 
maximum fidelity near 0.5. 

Fig’s[3ja) and (b) isolate the effect of exciton detuning 
and temperature, respectively, on the different fidelity 
characteristics. Fig’s [3](c) and (d) show the dephasing 
rate r f x over the same ranges, along with each of the 
contributing transition rates and the effective indirect 
exciton population decay rate T e //. As in Fig. [2j re¬ 
gions of high fidelity correlate with a radiatively-limited 
dephasing rate. The significant drop in fidelity at low 
exciton detunings is therefore due to a phonon-assisted 
hole tunneling rate which surpasses radiative emission by 
up to 3 orders of magnitude. Note that the oscillations in 
the tunneling rates are due to the phase appearing 
in the phonon coupling (Eq- ©), and depend on the 
distance d between QDsP2 The peak in fidelity width as 
a function of temperature occurs when T e jf increases, 
while "f d X remains low. The increased decay rate of the 
indirect exciton causes the time-integration in Eq. 0 
to sample shorter time delays, maintaining high fidelity 
over a larger range of exchange splittings. While hole 
tunneling can be adequately suppressed by simply main¬ 
taining a large enough detuning between exciton states 
or increasing the tunnel barrier, phonon-assisted spin-flip 
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FIG. 3: (a) (b) Dependence of fidelity parameters F^ ax 

(black, left y-axis) and A F + (red, right y-axis) on (a) de¬ 
tuning and (b) temperature, with th = 0.2 /reV. (c) (d) Total 
dephasing rate ~(' d of the indirect exciton bright state coher¬ 
ence, including contributions from incoherent transitions, as 
a function of detuning and temperature, respectively. Also 
included is the effective |zX) population decay rate r e //. 
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FIG. 4: (a) Maximum fidelity and (b) width of F + (Sd ) curve 
as a function of hole tunnel coupling and temperature, with 
A = 7 meV. 


surpasses radiative emission and reduces fidelity at tem¬ 
peratures above ~ 50 K. 

Fig- 0 shows the fidelity characteristics as a function 
of hole tunnel coupling and temperature, at a fixed ex¬ 
citon detuning. The effective indirect exciton decay rate 
T ef f depends strongly on both temperature and tun¬ 
nel coupling, and is demonstrated to have a substan¬ 
tial effect on the entanglement fidelity. F+ aa , remains 
higher than 0.9 where T e ff < 30Fjx, which occurs at 
low temperature and tunnel coupling values. The fi¬ 
delity width AFg peaks at intermediate values of T e // 
(T e ff ss 5 — lOlAx) and increases with decreasing tunnel 
coupling, reaching values above 3.0 peV at small tunnel 
couplings (th < -1 meV) and temperatures up to 70 K. 
As tunnel coupling increases, the indirect exciton gains 
more of a direct exciton character, decreasing the fidelity 
width and requiring lower temperatures to suppress re¬ 
combination. 


Our results indicate that by maintaining a large 
enough exciton detuning to suppress phonon-assisted 
tunneling, CQD-based entangled photon sources can pro¬ 
duce entangled photon pairs with higher fidelity than sin¬ 
gle QD-based sources, over a wider range of direct exciton 
exchange splittings. While this helps reduce the strict 
symmetry requirement for entanglement generation, the 
photons are separated by a relatively long time delay of 
100 — 1000 ns due to the low indirect exciton recombina¬ 
tion rate. In the intermediate state, however, the indirect 
exciton spin is entangled with the polarization of the first 
emitted photon (Eq. [2]). As a result, this scheme could 
potentially be used to entangle spins in remote CQDs by 
tuning each of the first emitted photons into resonance 
via electric or strain fields and performi ng a joint polar¬ 
ization correlation measurement; 17 EiHAl 


IV. CONCLUSION 

We have used a simple theoretical model to simulate 
the radiative cascade of the molecular biexciton state in 
a vertically stacked tunnel-coupled quantum dot pair. 
The entanglement fidelity of the resulting two-photon 
polarization state is determined, accounting for phonon- 
assisted tunneling and spin-flip processes. From numeri¬ 
cal simulations, we find an approximately Lorentzian de¬ 
pendence of fidelity on anisotropic electron-hole exchange 
splitting and mapped the behavior over a range of elec¬ 
tric field and temperature values. Our results show that 
near-unity maximum fidelity can be achieved over a range 
of exchange splittings \Sn\ <2.7 ^eV at large exciton de¬ 
tunings and low temperatures, where dephasing due to 
phonon-assisted hole tunneling and spin-flip processes is 
suppressed. This suggests that coupled quantum dots 
can generate photon pairs with a high degree of entan¬ 
glement and over a wider range of exchange splittings 
compared to single dots, provided the tunnel coupling is 
low enough to maintain charge separation in the indirect 
exciton state. In addition, the spin-photon entanglement 
generated by the first recombination could be used to 
entangle spins in remote CQDs. 
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Appendix: Computational Methods 

Electron and hole single-particle wavefunctions are cal¬ 
culated using the two-band effective mass approximation. 
We model the confinement potential of each QD as a fi¬ 
nite well in the growth (z) direction and parabolic con- 
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finement in the lateral (x-y) directions, expressed for par¬ 
ticle a = e,h as V a {r) = \m. a j n \ s io 2 a (x 2 + y 2 )+£ a 0(\z\ - 
|), where r7i a> i n As is the effective mass of particle a in 
InAs, ui a is the angular frequency of the lateral harmonic 
oscillator, and £ a is the conduction/valence band offset 
of the heterostructure. Envelope wavefunctions of single¬ 
particle localized ground states are then defined on a 
Cartesian grid as 


ip a (r) =Ae 


x 


(z 2 +y 2 ) x 

cos(kh/2)e~ Kiz - h/2) 
cos (kz) 

cos(kh/2)e K ( z+h l 2 ' > 


\iz> h/2 ( 14 ) 

if \z\ < h/2 
if z < — h/2 , 


where the wavenumber k is determined by the transcen¬ 
dental equation tan(fch/2) = 


/™^aT 7 m 77 ) 

W m Q ,GaAs J 


— 11 due to 


the boundary conditions, k = k\J — 
yj2m a j n A s £a/h, and A is a normalization constant de¬ 
fined such that / |VVb/t(0| 2 d 3 f = 1, and h is the QD 
height in the growth direction. In a coordinate system 
with the origin set at the center of the two QDs, wave- 
functions for particles localized in each dot are found 
from Eq. ([bl]) as ip a , B / T {x,y,z) = ip a (x,y, z±d/2), with 
the substitution h h B / T . The energies of these single¬ 
particle ground states are then given by 


h 2 k 2 

E a ,B/T = a- \- Huj a . (15) 

2 TO Qj InAs 


The Coulomb interaction energy between particles a 
and /3 in dots i and j, respectively, is given by 


y<*P _ J 3 _ 


d 3 rd 3< 


t, \^ai{r)\ 2 \^0oif')\ 2 


(16) 


where e is the charge of an electron, e is the average per¬ 
mittivity of InAs and GaAs, and the 6-dimensional inte¬ 
gral is directly evalutated by numerical integration over 
a Cartesian grid. The Coulomb interaction terms mod¬ 
ify the energy of multiparticle exciton states, as does the 
local electric field F due to the large permanent dipole 
moment p = ed of indirect exciton states. The energies 
of the various charge states considered are then 

Ex =E g + E eB + E hB — V BB 
Eix = E g + E eB + EhT — V B t — edF 
Eixx =2E g + 2 E eB + E hB + E hT 

+ V% e B + V$ - 2V^ h B - 2V& - edF, 


where E g is the band gap of the strained InAs comprising 
each QD. 

Eq. § is expanded using the interband dipole mo¬ 
ments = /it(±l, — i,0)/\/2 to obtain the expres- 


TABLE II: List of material and CQD parameters used in nu¬ 
merical simulations. 

Material Parameters 
Permittivity 
InAs electron mass 
GaAs electron mass 
InAs hole mass 
GaAs hole mass 
Strained InAs band gap 
Conduction band offset 
Valence band offset 
Hole deformation potential 
Mass density 
Speed of sound 
CQD Parameters 


Bottom QD height 


3.0 nm 

Top QD height 

hx 

2.5 nm 

QD separation 

d 

7.0 nm 

Electron level spacing 

hide 

66 meV 

Hole level spacing 

hw h 

15 meV 

Hole tunnel coupling 

th 

200 /mV 

Bright-dark splitting 

Am 

300 /mV 

Direct exciton recombination rate 

r x 

1.0 ns -1 


e 

13.5 e 0 

^e,InAs 

0.059 mo 

^e,GaAs 

0.042 mo 

^/i,InAs 

0.37 m 0 

^■/i,GaAs 

0.34 m 0 

e 3 

866 meV 


461 meV 

£h 

192 meV 

av 

700 meV 

p 

5300 kg/m' 

CLA 

5150 m/s 


Sion 


S D (I) = 


2ne 


d 3 fd 3 \ 


2Az 2 — Aa: 2 — Ay 2 + QiAxAy 
A r 5 

x ^lB{^eB{f')^l B(T) {r')4\ B(T ){r), 

(18) 

where Ar = r — r'. Since the exchange splitting varies 
widely between QDs and small values are required for 
entanglement generation, we choose values of S B and nu¬ 
merically integrate Eq. (18) to obtain the ratio Sj/S B - 
We use the value of Tx listed in Table [TT] and the 
electron-hole wavefunction overlap 

M ij = J d 3 rip* ei {r)ij> hj {r) (19) 


obtained by numerical integration to calculate T^y = 
F x(M B t/M bb ) 2 . The phonon-assisted tunneling rates 
I 'a its and T em aie calculated using Eq s |^bl| and 
|12[ where the phonon wave vector q is expressed in 
spherical coordinates and its magnitude is constrained 
by the delta function to be q = A/he la- To calcu¬ 
late the phonon-assisted spin-flip rate of the direct ex¬ 
citon, we use the experimentally-determined tempera¬ 
ture dependence T^ = (.27 ns -1 ) + (.29 ^s _1 K _2 )T 2 , 
and from it determine the indirect exciton spin-flip rate 
T l J = T^{M B t/M BB ) 2 . For the pure dephasing rate, 
we use the empirical temperature dependence fry* = 
(71 peV)n B (6 meV,T) + (4.5 meV)ns(28 meV,T) cor¬ 
responding to phonon-assisted pure dephasing processes 




















involving excited hole states, and assume the same pure 
dephasing rate for indirect excitons: 7 * x = 7 ^. 

By projecting Eq. [5] onto each combination of basis 
states, it becomes pij = Ylu ^and we determine 
the elements of the time-dependence tensor M tJ j,;. We 
then vectorize p by mapping the matrix elements pij to 
a one-dimensional vector p. The time-dependence ten¬ 
sor is then transformed into a matrix A'I 2 j , and 

the solution is expressed by a matrix exponential as 


p(r) = e Cr p( 0) = e Mt p( 0). With this matrix exponential 
describing time evolution of the CQD density matrix, we 
use Eq’s [ 6 ] and [7] to calculate the various elements of the 
two-photon polarization density matrix p 2 p in the lin¬ 
ear polarization basis {\HH), \HV ), \ VH ), |UU)}, aver¬ 
aging over delay times up to Td = 200 ns. With these 
density matrix elements, the fidelity to the entangled Bell 
state is finally calculated using Eq. [4j 
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